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Abstract. Acoustic Bessel beams are known to produce an axial radiation force on a 
sphere centered on the beam axis (on-axial configuration) that exhibits both "pressor" 
and "tractor" behaviors. The pressor and the tractor forces are oriented along the 
beam's direction of propagation and opposite to it, respectively. The behavior of the 
acoustic radiation force generated by Bessel beams when the sphere lies outside the 
beam's axis (off-axial configuration) is unknown. Using the 3D radiation force formulas 
given in terms of the partial wave expansion coefficients for the incident and scattered 
waves, both axial and transverse components of the force exerted on a silicone-oil sphere 
are obtained for a zero- and a first-order Bessel vortex beam. As the sphere departs 
from the beam's axis, the tractor force becomes weaker. Moreover, the behavior of 
the transverse radiation force field may vary with the sphere's size factor ka (where k 
is the wavenumber and a is the sphere radius). Both stable and unstable equilibrium 
regions around the beam's axis are found depending on ka values. These results are 
particularly important for the design of acoustical tractor beam devices operating with 
Bessel beams. 
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1. Introduction 

The transfer of momentum from electromagnetic [Ij and acoustical [i2j propagating waves 
upon reflection or absorption from a particle induces forces and torques that may be used 
under various circumstances to accelerate [3j, rotate ^ El El [7], trap [HIE], levitate [lOj 
or even stretch [11] the particle itself. The idea of attracting/repelling objects from a 
distance with radiation (known as tractor/pressor beams) has also captured the interest 
of many scientists and news media jT2l IT3l IT41 [T5j . The pressor behavior is somehow 
intuitive, such that the particle facing the incident beam is expected to be pushed 
away from the source. Alternatively, the tractor behavior occurs as a consequence of a 
negative force that pulls the particle toward the source. 

In acoustics, various investigations revealed effects similar to both pressor and 
tractor beams. For the pressor behavior, a large body of experimental and theoretical 
work has been already published (see Ref. |I6]). On the other hand, it was demonstrated 
that a curved- wave front, such as a progressive spherical diverging wave incident upon a 
rigid or solid sphere, induced a negative (attractive or tractor) axial radiation force 
depending on both the frequency of the incident fleld and the distance from the 
source jTTl [HI [IHl [20] . The same has been demonstrated for rigid cylindrical structures 
in a progressive cylindrical diverging waves [211 122]. A plane progressive wave can 
induce an axial negative force on an elastic sphere coated with a viscoelastic (absorptive) 
layer [23l [24]. Also, the presence of a boundary close to a spherical shell has shown to 
induce a negative axial acoustic radiation force that is aflFected by the frequency of 
the incident plane waves, the distance to the boundary and its corresponding reflection 
coeflicient [25j. In addition, when the incident beam is in the form of a zeroth-order 
(non- vortex) Bessel beam, the axial force can pull fluid spheres toward the source [26j. 
Moreover, the axial negative force of higher-order Bessel (vortex) beams on rigid |27j, 
fluid and elastic spheres [28j occurs depending on the beam's half-cone angle ^ and 
the sphere's size factor ka. An attractive radiation force can be also achieved by using 
highly focused acoustical tweezers p9j . 

In this work, the 3D radiation force caused by a zero- and a flrst-order Bessel 
beams upon an absorbing fluid sphere (silicone) is computed. In doing so, the analytical 
expressions for the axial and transverse components of the radiation force [30] are used. 
These expressions are given in terms of the beam-shape coeflicients of the incident wave 
and the scattering coeflicients of the sphere. The beam-shape coeflicients are evaluated 
through the discrete spherical harmonic transform [31 J developed in the context of 
the arbitrary acoustic scattering theory of Bessel beams by a sphere [32l EH]. The 
computations of the radiation force components for the silicone-oil sphere reveal that 
the attractive force of the tractor beam decays as the sphere departs from the beam's 
axis. In addition, the transverse radiation force fleld may exhibit stable and unstable 
equilibrium regions in the vicinity of the beam's axis. This study may potentially assist 
in designing speciflc acoustical probes operating with Bessel beams for biological imaging 
and particle manipulation applications. 
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2. Acoustic radiation force 



Consider an incident acoustic wave of angular frequency u scattered by a fluid sphere of 
radius a suspended in the wavepath within an ideal fluid. The pressure of the incident 



wave is given by pi 



Pie 



-iujt 



, where pi is the pressure amplitude, i is the imaginary 



unit, and t is the time. In turn, the fluid is characterized by the ambient density po 
and the adiabatic speed of sound Cq The 3D radiation force generated on a suspended 
sphere can be written in the form f = 7ra^£'o(Y^, 1^, Y^), where is the characteristic 
energy density of the wave, Y^^ Yy and are the radiation force functions in Cartesian 
coordinates. These functions are given by [30\ 

= Re{Y^}, (1) 

Yy = lm{Y^}, (2) 
1 



Y^ = 



7r{ka)' 



(3) 



Lm 



where and s]^ are the beam-shape and the scattering coeflicients, the symbol * means 
complex conjugation, = \/ {I — m){l + m)/[{2l — 1)(2/ + 1)], 'Re' and 'Im' denote the 
real and the imaginary parts, and 



Yi 



^ Lm 



^)(c^rn-li-m 



(4) 



with b]^ = ^{l + m){l + m + l)/[(2/ - 1)(2/ + 1)]. 

The beam-shape coeflicients are determined from the incident pressure amplitude 
as follows [32l [33] 

1 



p,{kb,9,^)Yr{9^^)dn, 



(5) 



Poji{kb) 

where po is the incident pressure magnitude, ji is the spherical Bessel function of l-th 
order and dQ = sin 9 d9 dip. The integral is performed over the surface of a control 
spherical region of radius b. This region encloses the inclusion that is subjected to the 
radiation force. The integral in Eq. ([5]) is the spherical harmonic transform (SHT) of 
the incident pressure amplitude. For an arbitrary pressure amplitude the integral in 
equation ([5j) should be evaluated numerically. 

The scattering coeflicient for an absorbing fluid sphere can obtained from the 
following boundary conditions: the pressure and the radial component of the particle 
velocity should be continuous across the sphere's surface. Accordingly, we have [7J 

-fji{ka) ji{kia) --fh[^\ka) ji{kia) 

j[{ka) j[{k^a) \ ^ [ -hf^\ka) j[{kra) 

where 7 = kipo/kpi is the impedance index, ki = uj/ci + ia^ with a being the absorption 
coeflicient, /i^^ is the flrst-kind spherical Hankel of l-th order, and the prime symbol 
means diflFerentiation. 



a, 



det 



(6) 
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3. Acoustic Bessel beam 

The pressure amplitude of a n-th order Bessel beam propagating along the z-axis is 
expressed as ^34] 

Pi,n = e^'^^°^^Jn(fc^sin/5)e-^ (7) 

where is the nth-order Bessel function, /3 is the half-cone angle of the beam, 
Q = \^{x — XqY + {y — y^Y-) and Lp = tan~-^[(^ — yo)/{x — Xq)]. The relative position of 
the beam axis with respect to the center of the sphere in the x^-plane is denoted by 
(^0,2/0). 

The beam-shape coefficients for an off-axial Bessel beam need to be numerically 
computed from the SHT given in equation (Ej). This is performed through the discrete 
SHT (DSHT), which is realized through the fast Fourier transform (FFT) for the 
integration in the polar angle and the Gauss-Legendre quadrature used for the azimuthal 
angle integration [31j. In this algorithm, the pressure amplitude is sampled on the 
control sphere surface. In the azimuthal direction, this function has M uniformly 
distributed sampling points, while in the polar direction, it has points corresponding 
to the zeros of the Legendre polynomial of A^-th order. 

4. Results and discussion 

The radiation force caused by a zero- and a first-order Bessel beams upon a sphere 
suspended in water (po = 1000 kg/m^ and Cq = 1500 m/s) is analyzed here. The 
sphere is made of silicone oil [35j (pi = 970kg/m^, Ci = 1004 m/s, and a = 
2 X lO"-*^^/-^" Np • Hz~-'^ ''/m, with / being the frequency). The silicone oil was chosen 
because it may form drops into water without mixing. The frequency of the incident 
beams is fixed to 1 MHz. 

To compute the radiation force, we have to obtain the beam-shape and then the 
scattering coefficients given in equation ([6]). To compute the beam-shape coefficient we 
use the DSHT algorithm with M = 2^^ and = 256. In this case, the mean absolute 
error between the pressure amplitude of the Bessel beams given in equation ([71) and 
reconstructed waves using the beam-shape coefficients is about 10~^. Furthermore, 
the infinite series in equations ([3]) and (|4j) should be truncated at / = L. The 
truncation order L should be the smallest integer for what the scattering coefficient 
ratio I^^/sqI <C 1. This is required to ensure proper convergence of the radiation force 
series. In the performed computations, the truncation order is such that the scattering 
coefficient ratio is as small as 10~^. 

To analyze the behavior of the axial radiation force, 2D plots of the radiation force 
Yz are computed in the range 0° < /3 < 90° and < fca < 10. As an initial test, we 
have recovered previous results based on the analytical calculation of the axial radiation 
force upon an hexane sphere in the on- axis configuration [28]. For the sake of brevity, 
these results will not be shown here. 
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The 2D plots of the axial radiation force function produced by the zero-order 
Bessel beam are depicted in figure [H The radiation force function is normalized to 



1.3931. The on- axis configuration is shown in figures 1(a) and |l(b) , In this configuration, 
islands of negative radiation force (tractor beam) are clearly spotted. The largest 
magnitude of the negative axial force lies in the region around ka = 2 and 60° < /3 < 80°. 
In this region, the attractive force magnitude corresponds to approximately 8% of the 
maximum value of the repulsive axial force. With an oflFset of kxo = 1.6 along the x-axis, 
the islands of attractive force shift toward higher ka values. Moreover, the strength of 
the negative force drops one order of magnitude. 
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Figure 1. Axial radiation force function for the silicone oil sphere due to the zero- 
order Bessel beam. The function Yz is normalized to 1.3931. (a) The sphere is placed 
on the axis of the beam, (b) The islands of attractive force in the on-axis configuration, 
(c) The sphere location has a offset of kxo = 1.6 along the x-axis. (d) The islands of 
attractive force in the off-axial configuration. 



Results of the axial radiation force function due to the first-order Bessel vortex 
beam as a function of /3 and ka are shown in figure [21 The radiation force function 
is normalized to 0.4201. In the on-axis configuration (figures [2(a) and 2(b)D , islands 
of attractive radiation force arise. The largest magnitude of the negative force takes 
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place within the region 65° < /3 < 80° and ka = 3. The attractive force here is about 
20% lower than that due to the zero-order Bessel beam. When an offset of kxo = 1.6 
is applied to beam along the x-axis, the islands of attractive force move toward higher 
ka values. Furthermore, with this offset the attractive force islands shrink and their 
strength drop by one order of magnitude. 




(c) (d) 

Figure 2. Axial radiation force function Yz for the silicon oil sphere due to the 
first-order Bessel beam. The function Yz is normalized to 0.4201. (a) The sphere is 
placed on the beam's axis, (b) The islands of attractive radiation force in the on-axis 
configuration, (c) The sphere location has a offset of kxo = 1.6 along the x-axis. (d) 
The islands of attractive radiation force in the off-axial configuration. 

The transverse radiation force field = (Yx^Yy) caused by the zero-order Bessel 
beam with /3 = 70° for the silicone sphere is illustrated in figure [31 The transverse 
force fields are computed by varying the relative position between the sphere and the 
beam (offset) within —6.4 < kx^^kyo < 6.4. Furthermore, the corresponding vector 
fields are plotted on top of the time-averaged energy flux (intensity) of the incident 
beam normalized to the unit. Two different sphere's size factors are considered, which 
correspond to a repulsive {ka = 0.1) and an attractive {ka = 2) axial radiation force (see 
figure [I]). In both cases, the central spot of the beam is a region of unstable equilibrium 
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for the sphere. The hnes of force diverge from the beam axis. 

In figure [H the transverse radiation force field due to the first-order Bessel beam 
with (5 = 70° for the silicone sphere is depicted. The force fields are computed similarly 
to those in figure [3l To analyze the pressor and the tractor beam two size factors are 
chosen, namely fca = 0.1 and ka = 3. For fca = 0.1 the beam central spot is a region of 
stable equilibrium for the sphere. The sphere is pushed by the pressor beam along the 
z-axis. The transverse force points inwardly toward the beam axis. Thus, the sphere 
can be trapped in the beam's axis. As the sphere is moved radially away from the 
beam's axis, the transverse radiation force reverses its direction. It fiips from inward to 
outward direction as the spheres varies from the first maximum to the first minimum 
of the beam's intensity. Moreover, the transverse radiation force points inwardly in the 
region from the first minimum to the second maximum beam's intensity. In this region, 
the sphere might be trapped. The transverse radiation force field changes considerably 
when ka = 3. It points outwardly following the clockwise direction as the sphere moves 
radially away from the beam's axis. Moreover, inside the beam's central spot, the sphere 
is pulled back along the z-axis by the tractor beam. Moving the sphere radially from 
the first maximum to the first minimum, the transverse radiation force points outwardly 
along the counterclockwise direction. No stable equilibrium region is formed in this case. 



5. Conclusion 

The radiation force produced by a zero- and a first-order Bessel beam upon a silicone- 
oil sphere in both on- and ofl[-axial configurations has been studied. The analysis was 
based on the 3D radiation force expressions given in terms of the beam-shape and the 
scattering coefficients. The beam-shape coefficients were numerically computed through 
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Figure 4. Transverse radiation force field = (Yx^Yy) caused by the first-order 
Bessel beam with /3 = 70° upon the silicone sphere with (a) ka = 0.1 and (b) ka = 3. 
The background images correspond to the beam intensity. 



the DSHT algorithm, while the scattering coefficients were obtained by applying the 
acoustic boundary conditions across the sphere's surface. With appropriate selection of 
the sphere size factor ka and the beam half-cone angle /3, the axial radiation force exerted 
on the sphere might be repulsive (pressor beam) or attractive (tractor beam). As the 
sphere departs radially from the beam's axis, the attractive axial radiation force becomes 
weaker. In addition, the transverse radiation force field was computed showing that 
the silicone sphere can be trapped by using the first-order Bessel beam. These results 
represent an important step toward the design of a new generation of acoustical tweezers 
operating with Bessel beams for potential applications in bioengineering, biophysics, and 
other related fields. 
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